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Abstract 

As soon as the geometry expands quasi-exponentially the plasma sources are customarily 
tuned to zero as if the duration of the inflationary phase was immaterial for the gauge 
field fluctuations at large-scale. The serendipitous disappearance of the plasma (or even 
the partial neglect of its effects) depends on the symmetries of the system which are, in 
four space-time dimensions and in the simplest Abelian case, the invariance under Weyl 
rescaling and the electromagnetic duality symmetry. The quantum, thermal and conducting 
initial conditions of inflationary magnetogenesis are classified and discussed with the aim of 
determining when plasma effects can be effectively disregarded. The speculative implications 
of a non-degenerate monopole plasma for the conservation of the large-scale electric flux are 
briefly examined. 
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1 Introduction 



In four-dimensional curved space-times, besides the invariance under local gauge transforma- 
tions, two symmetries are directly relevant to the gauge fields: the Weyl symmetry jl] and 
the electromagnetic duality symmetry [21 13] . If the governing equations of a given quantum 
field are invariant under Weyl rescaling, the corresponding normal modes are not excited by 
the evolution of the geometry [U |5l [6] and particles are not produced. Duality rotates field 
strengths into their duals (i.e. tensors into pseudotensors) and it is therefore not an internal 
symmetry. 

The symmetries of the evolution equations of the gauge fields in curved space-times are 
the main handle on the origin of the large-scale magnetism [3 |H1 E], a problem dubbed 
magnetogenesis some time ago [10] but whose origin is much older as a number of relatively 
ancient but still very inspiring monographs demonstrates [HI [121 US]- Iii the traditional 
framework magnetogenesis can take place either during or after a conventional phase of 
inflationary expansion. For this reason inflationary and post-inflationary magnetogenesis 
have been regarded as two separate and mutually exclusive possibilities. Amplified quantum 
mechanical fluctuations are associated with inflationary magnetogenesis by implicitly dis- 
regarding the role of electromagnetic sources. Conversely, post-inflationary magnetogenesis 
rests on the existence of appropriate currents operating inside the particle horizon (if present) 
at a given epoch in the life of the Universe. The purpose of this paper is to bridge the two 
aforementioned perspectives by scrutinizing the role of plasma sources in the situation where 
they are customarily neglected, i.e. during inflation. As it will be shown, the symmetries 
of inflationary magnetogenesis do not forbid the presence of electromagnetic sources unless 
they are fine-tuned to vanish initially for some theoretical prejudice. The present analysis 
completes and systematizes the results reported in [13] where a swift account of some of the 
ideas contained in this paper has been illustrated. 

The duration of inflation is parametrized in terms of the total number of efold^, denoted 
by A^. Since A^ is not accessible to direct observations it is customary to introduce various 
(conceptually related) quantities such as A^min (i-e. the minimal number of efolds necessary 
to fix the problems of the standard big-bang cosmology) or Amax (i-e. the maximum value 
of inflationary efolds presently accessible to our cosmological observations). The value of 
Amin stems directly from the analysis of all the kinematical and dynamical problems of the 
standard hot big-bang scenario while the value of A^max can be determined by fitting the 
present size of the Hubble radius inside the event horizon of the quasi-de Sitter space-time. 
The numerical value of Amax is not independent on the post-inflationary thermal history 
and may range from about 62 (for the standard thermal history with sudden reheating) 
to larger figures if the reheating is delayed (see, for instance, [151 II^)- For the standard 

^The number of efolds is the natural logarithm of the ratio between the scale factor at end (i.e. Qc) and 
the beginning (i.e. Ob) of the inflationary expansion. For the purposes of this introductory section we can 
conventionally agree that inflation starts as soon as an event horizon is formed. 
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thermal history with sudden reheating A^min — ^max- In the case of scalar and tensor modes 
of the geometry, the conventional argument stipulates that when the total number of efolds 

is much larger than iVmax the effect of the initial conditions matters very little for any 
conclusion involving the amplification of quantum fluctuations. In the limit ^ iVmax the 
quantum fluctuations are automatically identified, in practice, with vacuum fluctuations. 
The problem of the initial conditions for the scalar and tensor fluctuations of the geometry 
has been discussed in the past; see, for instance, Refs. [T71 [HI |T9] and [201 l21l |22l |23] for a 
still incomplete list of articles. 

In contrast to what happens for the scalar and tensor modes of the geometry, it will be 
argued that the Weyl and duality symmetries determine if and when the initial conditions 
of the problem are to be treated classically or quantum mechanically. Even conceding that 
the total duration of inflation is unknown, the inflationary phase is not eternal in the past 
either and potential sources cannot be neglected at the protoinflationary transition, i.e. 
when the background geometry starts accelerating. An exceedingly large number of efolds 
does not justify, per se, the total or partial neglect of the electromagnetic sources during 
the inflationary evolution of the gauge fields. If the protoinflationary initial conditions 
are dominated by a relativistic plasma (containing either electric or magnetic sources) the 
whole system is invariant under Weyl rescaling. In the latter case the conductivity, may 
stay unsuppressed up to a critical efold (be it A^c) which is not determined by the initial 
conditions but by the mass of the charge carriers. Charge carriers sufficiently light during 
inflation and with masses in the GeV range lead to A'c ~ C(34) (see, for instance, section H]). 
In this example the sources can still be relevant (A^max — A^c) efolds before the end of inflation 
in contrast to the conventional set of assumptions. It is therefore not excluded that Weyl 
symmetry prevents the dissipation of protoinflationary Ohmic currents by so modifying the 
standard protocol for the assignment of the initial data in inflationary magnetogenesis. 

In a realistic approach to the quantum and classical initial conditions of inflationary 
magnetogenesis it is safe to enforce the global charge neutrality of the primeval plasma prior 
to the onset of inflation to avoid any phenomenological drawback [23]. This requirement 
implies that the concentration of positively charged species is exactly balanced by the neg- 
atively charged ones. Denoting with ?t,+ and n_ the comoving concentrations of the species 
with positive and negative electric charges, we shall assume that n+ = n_ = no where no 
is the comoving value of the common charge concentration determining directly the Debye 
shielding length of the initial electric field. If no 7^ during the protoinflationary phase the 
evolution of the electric and of the magnetic fields is sharply different. The presence of an 
electric current breaks the duality symmetry which is instead unbroken in the absence of 
electromagnetic sources. The latter situation leads to the classical (or conducting) initial 
conditions. Conversely, if the charge concentration vanishes, i.e. no = 0, two complemen- 
tary subclasses of initial conditions may arise depending on the duration of inflation. If the 
duration of inflation greatly exceeds A^max the initial conditions minimize, in practice, the 
Hamiltonian of the gauge field fluctuations. Finally if the total number of efolds is close to 
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-^max (but always in the case Uq = 0) the initial conditions for the evolution of the electric 
and magnetic fields are likely to be thermal. 

The present paper is organized as follows. Section [2] introduces the Weyl and duality 
symmetries in conformally flat background geometries. In section [3] the quantum mechanical 
treatment of the vacuum and of the thermal initial conditions is addressed when the total 
charge concentration vanishes (i.e. uq = 0). In section H] the plasma initial conditions, 
corresponding to the case Uq ^ 0, are examined. Section [5] contains some considerations on a 
Lorentzian plasma of magnetic monopoles. The concluding remarks are collected in section 
El To avoid lengthy digressions, some among the most relevant technical results have been 
collected in the appendix. 



2 Weyl and duality symmetries 

The standard procedure for assigning the initial data to large-scale fluctuations stipulates, 
in a nutshell, that as soon as the background geometry inflates, all the sources present in the 
plasma can be serendipitously neglected with the sole exception of an (effective) cosmological 
term mainly produced by the potential of the inflaton. The previous recipe is justifled when 
applied to the evolution of the geometry itself but it is arbitrary as far as the governing 
equations of the gauge flelds are concerned. The point is that given a system of equations, 
there can be symmetries preventing the dissipation of the corresponding sources. To shed 
light on this aspect, the Weyl symmetry and the duality symmetry will now be swiftly 
scrutinized. 



2.1 Weyl symmetry 

Consider the following rescaling of the four- dimensional metric tensoilfl 

G^,u gf,u = q{x) G^u, x = x'' = (r, x'), (2.1) 

where x denotes the space-time coordinate. Under the transformation (12.1 1) a generic Abelian 
vector field, the related field strength and its dual become: 

y^^^Y^ = 4-, 3^^^ ^ Y^^'' = 4—, 3^^^ ^ r^'^ = 4—, (2.2) 
q[x) q^[x) 

where Yap and Y'^^ are, respectively, the field strength and its four- dimensional dual defined 
with respect to the metric g^^. 

Yo.fi = V[„y^] = d[aYp], Y-" = -^Yp,- (2.3) 



-9 



•^The signature of the metric is taken to be mostly minus. Four-dimensional indices are denoted by 
lowercase Greek letters while the spatial indices are denoted with lowercase Latin characters. Units h = c = 
Kb = 1 will be assumed hereunder. 
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Vq, is the covariant derivative with respect to the metric g^^; yap and 3^"^ are the field 
strength and its dual but defined in terms of the metric G^u- Using Eqs. (12. ip and fl2.2p it 
follows, in four space-time dimensions, that the vector potential transforms as: 



V^y" V^Yf" = q{x) ^^y^", (2.4) 
while the corresponding field strengths transform as: 



v=G y^"" V^Yf"" = 3^^^ y"" ^ y^^ = 3^^^ (2.5) 

Equation (12. 5p implies that \/—g Y^^ and ^/^^Y^^ are both invariant under Weyl rescaling. 
The generally covariant four-divergence of the field strength and of its dual vanish in the 
vacuum (i.e. Y^^'' = Y'^'^ = 0) as a consequence of the Maxwell's equations that can 
also be expressed as d^{^—g Y^^) = d^{y/—g Y^^) = 0, provided det( ^f^i, ) 7^ 0. The previous 
chain of arguments proves the invariance of the original equations under Weyl rescaling 
thanks to Eq. (12. 5p . The generally covariant Lorentz gauge condition (i.e. VpY^ = 0) is 
not Weyl invariant as it follows from Eq. (12. 4p . Consequently, the Coulomb gauge turns out 
to be more suitable in curved space-times. This observation will be used in section |3] and in 
appendix \M 

Consider now, for sake of simplicity, the physical situation where G^i, = rj^y {rj^^ = 
diag(l, —1, —1, —1) denotes the Minkowski metric). Such a situation is qualified as physical 
since, whenever q{x) does not depend on the spatial coordinates but only upon the conformal 
time coordinate r, g^^^, describes, according to Eq. (12. ip . a Friedmann- Robert son- Walker 
geometry with fiat spatial sections. The simultaneous presence of gravitating electric and 
magnetic sources together with an effective cosmological term implies that the dynamics of 
gfj^u is dictated by the Einstein equations appropriately written in their contracted form: 



(2.6) 



where R^i, is the Ricci tensor computed from the metric g^y while T^y denotes the correspond- 
ing (total) energy-momentum tensor given as the sum of the energy momentum tensors of 
the sources (i.e. t^^), of the electromagnetic field (i.e. T^u) and of the (effective) cosmological 
constant A: ^ 

Tuu = tiiu + Tiiu + — 'p;giiu- (2-7) 
ovrG 

Since T^j^v is traceless, only the two remaining terms of Eq. (12. 7p will contribute to T = 
in Eq. (12. 6p . The Ricci tensor, appearing at the left hand side of Eq. (12. 6p reads, in explicit 
terms: 



R^u 2 



where Q = Inq. The scaling of Eq. (12.81) with Q must be compared with the energy- 
momentum tensor appearing at the right hand side of Eq. (12. 6p whose different contributions 
scale as follows: ^ 

oTTCt 
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In Eq. fl2.9p t^,^ and T denote the rescaled energy-momentum tensors of the sources and 
of the electromagnetic field. Since under Weyl rescaling both ^J—gY^^ and ^J—gY^" are 
invariant, must scale as = . Relativistic (or non-relativistic) charge carriers lead 
to different values of b in Eq. (12. 9p : quite generically, however, 5 > 0. In summary Eqs. (12.81) 
and (12. 9p show that Eq. (12. 6p is not Weyl-invariant. Furthermore the scaling properties of 
the sources at the right hand side of Eq. (12. 6p guarantee that as soon as Q increases, the 
effective cosmological constant dominates, at least asymptotically. 

The evolution equations of the gauge fields in the presence of electric and magnetic 
currents are, respectively 

V^F^^ = 47rj^ V^,Y^'' = Aixt, (2.10) 

we do not consider here dyon^ but either a plasma of electric charges or a plasma of magnetic 
charges. If = and ^ the plasma is purely electric since is a current of electric 
charges. If 7^ and j'^ = the plasma is purely magnetic since i'^ is a current of magnetic 
monopoles. As established from the general discussion of Eqs. (I2.4p - (l2.5p . Eqs. (I2.10p and 
(12. 7p are Weyl-invariant provided \/—gj^ or ^—gi^ are separately Weyl-invariant. Since we 
do not have any evidence of magnetic monopoles the electric plasma must be regarded as 
physically more realistic from the viewpoint of protoinflationary initial conditions. In the 
case of a purely electric plasma characterized by an Ohmic current Eq. (I2.10p can be 
written as: 



47rv^(T(a;)F^"u„, (2.11) 



together with the supplementary conditions dp,{\f—gY^^) = and g'^^ u^up = 1. Equation 
(12. lip is invariant under Weyl rescaling provided 

cr(x) — )• a{x) = Jq{x)(T{x), Ua{x) — )■ Mq,(x) = " (2-12) 

'q{x) 



The scaling law of the conductivity reported in Eq. (I2.12p holds, for instance, in the case of a 
relativistic plasma in approximate thermal equilibrium when the temperature exceeds both 
the mass of the charge carriers and the chemical potentials. This conclusion is well known [1] 
and has been used in a number of interesting analyses [231 [251 EZl I2B] but all related to the 
post-inflationary evolution of the electromagnetic fields. We notice here that the conclusions 
of Eqs. (12. lip and (I2.12p are qualitatively different from the ones deduced from Eqs. (12. 8p 
and (12. 9p : while asymptotically for large Q the electromagnetic sources are exponentially 
suppressed in the Einstein equations, they are not suppressed in Eq. (12. lip . Even if the 
electromagnetic sources can be ignored in the asymptotic dynamics of the geometry their 
presence cannot be ignored when assigning the initial data for the gauge fields. 



^Dyons are idealized point particles carrying both magnetic and electric charges; in this paper we shall 
confine our attention to the case of electric currents since this is, after all, the only realistic case. The 



magnetic currents are briefly discussed in section O 
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The symmetric argument can be repeated for a purely magnetic plasma (i.e. f'' 7^ and 
j'^ = 0) where Eq. f l2.10p becomes: 

= 47rv^(T^Z""n„, (2.13) 

together with the conditions df^{y/^ Z'^") = and g""^ UaUp = 1; in Eq. f l2.13p cr^ denotes 
the magnetic conductivity. As in the case of Eqs. (12.111) . y/—g Z^^ and ^—g Z^^ are sepa- 
rately invariant under Weyl rescaling. The right hand side of Eq. (I2.13p is Weyl invariant 
provided am and Ua scale, respectively, as \J q{x) and l/^Jq{x). A magnetic plasma is physi- 
cally very different from the case of an electric plasma. The passage from a purely electric to 
a purely magnetic plasma is an example of duality rotation which will be introduced hereun- 
der. The monopole plasma will be briefly discussed in section |5] since it is per se interesting 
but not strictly central to the main theme of this paper. 

The main lesson to be drawn is, in short, the following. The results of Eqs. (12. lip . 
(I2.12P and (I2.13P provide two examples where the Weyl symmetry prevents the damping of 
the electromagnetic sources. The sources cannot be neglected in spite of the duration of 
inflation unless they were already absent at the beginning or unless they are fine-tuned to 
zero by fiat. Even more concretely, if the temperature of the charge carriers is approximately 
comparable with the radiation temperature (as it happens when the plasma parameter is 
very small, see section H]) the whole plasma is Weyl invariant in the relativistic regime when 
the masses of the charge carriers and the chemical potentials can be neglected. 

It is appropriate to mention, as a side remark, that some of the considerations of the 
present section can be formulated in terms of an antisymmetric six-tensor and a six-current 
as suggested originally by Dirac [29]. The latter approach does not add any immediate 
advantage for the physical considerations of this analysis but has the virtue of enlightening 
the connections between Weyl spaces and the special conformal transformations [301 l3T] . 
These interesting connections will be left aside; we shall instead proceed with our four- 
dimensional formalism without grouping the field strengths and their duals in a (constrained) 
six- dimensional antisymmetric tensor field. 

2.2 Duality symmetry 

Consider a rotation of the field strength and of its dual parametrized by an angle "i?: 

Yf"" ^ Z^"" = cos'dY^"' + sm^Yf"', Yf"" ^ Z^"" = - sin^Yf"" + cos^Yf"" . (2.14) 
Under Eq. (I2.14p Eq. (I2.10p transforms as: 

V^Z^"" = AttL", V^Z^"" = AnJ", (2.15) 

provided 

f ^ r = 008^3" + sm^r, ^ L" = -sin^f + cos^r. (2.16) 
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The duahty transformation mapping an electric plasma into a magnetic plasma is realized for 
■(9 = 7r/2 when Z^"" = ^ = -F'^^ f ^ r = t and t = -f. 

In the latter case Eqs. fl2.12p and f l2.13p can be explicitly related. 

The duality symmetry can be generalized to the case when the kinetic term of the gauge 
fields is coupled to a scalar degree of freedom \{x). Consider then the following action: 

Sy = -T^ / d'xV^X{x)Y^pY''^. (2.17) 

Recalling that Y^ = (Yq, —Yi), in the Coulomb gauge (i.e. Yq = and V ■ F = 0), we shall 
have 

Yo, = a^e, = -drY,, F,, = -a^e.^^fo'^ = d.^^ - djY,. (2.18) 
The electric and the magnetic fields e and b are related to the canonical fields E and B as: 





E = d—a'e = -y' + J^y, B = \—aH = Vx y, (2.19) 



where J-" = / is the rate of variation of y/X. The consistent use of 5, E and y avoids 
potential problems in the identification of the electric and magnetic degrees of freedom. The 
canonical fields appearing in Eq. f l2.20p are exactly the normal modes of Eq. f l2.19p . In this 
source-free case the analog of Eq. f l2.10p can be written in terms of the canonical electric 
and magnetic fields: 

VM^XB)-I-XXE)^,, _^(^)Wx(|)^0. (2.20) 

The system fl2.20p is left invariant by the generalized duality transformation: 

E^-B, B^E, VA^-^. (2.21) 

The symmetry of Eq. ( ]2.2ip holds in the absence of electromagnetic sources. The separate 
addition of either an electric current or of a monopole current breaks the duality symmetry. 

In general terms A = \[ip{x), 4'{x), ...] may be a functional of various scalar degrees of 
freedom present in the model such as the infiaton ip [32], the dilaton [33], a dynamic gauge 
coupling [311 ES] (see also [361 EZ] cind [TO])- The field A can be a functional of a spectator 
field ip, [SSlEn] (see also [101 HI] ) evolving during the infiationary phase; in this case there is 
no connection between the evolution of A and the gauge coupling. The physical features of 
the various models are different: while some of these ideas are preferentially realized in the 
case of bouncing models [TOl |33] , some other are compatible with the standard infiationary 
paradigm [321 131 [351 138] . 

In the present investigation we shall preferentially consider models compatible with the 
conventional infiationary scenario where A depends on a spectator field. It is also possible 
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to think in terms of an effective gauge coupling but, in tliis case, it is necessary to assume 
that the gauge couphng is stronger at the beginning, then it decreases and so that y/X is 
bound to increase. The latter requirement fits perfectly with the observation that in the 
conventional inflationary expansion the gravitational coupling is strong: the curvature scale 
at which inflation takes place is 0(10"^ Mp) and, as we go backward in the past, the curvature 
gets even stronger. In this situation the electric fields are always deamplified and potential 
backreaction effects are negligible as discussed within different perspectives in |3H|35l|38]. If 
the gravitational coupling is weak at the beginning (i.e. the curvature is vanishing small and 
the initial conditions are given in flat-space time) it is natural to ask that the gauge coupling 
is also small at the beginning and this happens, for instance, in the case of bouncing models 
[SSI 1121113]. 

When the initial conditions are dominated by a globally neutral plasma the situation is 
yet different insofar as the correct expansion parameter which guarantees the validity of the 
plasma approximation is the plasma parameter (i.e. the inverse of the number of particles 
present in the Debye sphere) and not simply the gauge coupling constant. The plasma 
parameter is bound to decrease during the inflationary phase since this evolution guarantees 
that the Debye shielding scale gets larger, the electric flelds are screened and, in this way, 
they lead to irrelevant backreaction effects throughout the whole inflationary stage. 



3 Thermal and quantum initial condition 

In this section it will be assumed that the concentration of the positively and negatively 
charged species is flne tuned to vanish even if, according to the considerations on the Weyl 
and duality symmetries reported in section HI this is not the most generic situation. It 
is though interesting to see what happens when the conventional quasi-de Sitter phase is 
preceded by a decelerated stage of expansion 

lim ait) = lim d(t) > 0, lim ii(t) < 0, lim 'd(t) > 0, (3.1) 

where H^, deflnes the time-scale of the protoinflationary transition and a{t) is the scale 
factor in cosmic time. According to Eq. (13. ip . the acceleration changes its sign as soon 
as the inflationary event horizon is formed while the background expands both before and 
after the transition. The background geometry is conformally flat (i.e. Qf^u = 0''^{''~)i]^u) and 
falls into the more general class of examples analyzed in section [2j The connection between 
the cosmic time coordinate t and the conformal time coordinate r is given, as usual, by 
a(r) dr = dt. To avoid potential misunderstandings, it should be clear that the transition 
from a deceleration to acceleration has nothing to do with the so-called bouncing behaviour 
where the background passes from contraction to expansion or vice versa (see, e.g. |l2l US] 
and references therein). 
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3.1 Protoinflationary initial data 



In what follows, the protoinflationary sources are supposed to be in thermal equilibrium at 
a physical temperature such thato: 

HrMp = ^Pr^ Pr = ^^thTj, (3.2) 

where gth is the effective number of spin degrees of freedom. For if^, < Hr the universe 
inflates driven, for instance, by a single inflaton field so that the inflationary curvature scale 
is determined as 

where V denotes the inflaton potential and ^ is fixed by the amplitude of the power spectrum 
of curvature perturbations An at the pivot scale kp = 0.002 Mpc~^. In terms of the WMAP 
data alone [H] and, in particular, for the 7yr data release |15] analyzed in the light of the 
vanilla ACDM paradigrqj we have: 



^=^TTeAn, = (2.43 ±0.11) X 10"^ (3.4) 

where e = —H/H^ is the slow roll parameter. The present value of is 



TrM-Q.T,„,^[i^) (3.5) 
with _ 

Tr{tr) ^ f 45 



Note that < 1 as long as Pritr) < 3iJ^ Mp/ (Stt). The appropriate initial state is therefore 
a mixture characterized by a density matrix 

P = EAn}IM)(WI, An} = n .. , (3-7) 

{n} k ^ ~^ 

where A^^ is the average occupation number of each Fourier mode and, following the standard 
notation, |{n}) = ® ® Wk^--- where the ellipses stand for all the occupied modes 
of the field. 



^Temperature-dependent phase transitions [T71 [THl [HI HOI [211 HH 123] lead to an initial thermal state for 
the metric perturbations. The same kind of approach will be pursued in the present section. If the initial 
state is not thermal (but it is not the vacuum either) the present formalism applies with the difference that 
the spectral dependence of the mean number of photons will not necessarily have a Bose-Einstein form. 

^We remind that in the acronym ACDM the A stands for the late-time dark energy component and 
CDM stands for the late time cold dark matter component. The early completion of the ACDM paradigm 
contemplates a conventional inflationary phase with standard post-inflationary thermal hystory. 
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3.2 Two-point functions 

Recalling the results of appendix [Aj and, in particular, Eq. (lA.Sp . the canonical fields 
appearing in the duality-invariant Hamiltonian 



7r 



2J^Ti-y + diy-d'y 



(3. 



can be promoted to quantum operators (i.e. yi — j- iji and tTj — )■ tTj) obeying (equal time) 
commutation relations: 



[yi(xi,r),7rj(x2,r)] = i/\ij{xi - fa), 



Ai,-(fi - X2) 



.e'H^i-^2)p^^Q^^^ (3.9) 



where Pij{k) = {6ij — kikj/k"^). The function Ajj(xi — X2) is the transverse generalization of 
the Dirac delta function and such an extension is mandatory since V ■ -E = (because of the 
Gauss constraint) and V ■ y = (because of the gauge condition). The field operators will 
then become (see also Eqs. ( IA.4l) -( lAl5l) of appendix [A]) 



Tri[X,T 



(27r)3/2 ^ 



it e'''-^ 

k, a 



d^k 

(27r)3/2 



Y.^\k) 



9k{r)dre (7^(r) a 



it e'^-^ 

k, a 



(3.10) 
(3.11) 



where, e\°'\k) (with a = 1, 2) are two mutually orthogonal unit vectors which are also 
orthogonal to k; the creation and annihilation operators appearing in Eqs. f l3.10p and (13. lip 



obey [d 



^ais^^^Kk + p). Fom Eq. flA.7p the mode functions fk and obey: 



f'k = 9k + J^fk 



a'k = -k^fk - J^gk 



(3.12) 



and satisfy the Wronskian normalization condition fk{'^)gli'^) ~ fk{'^)9k{'^) ~ which follows 
from enforcing the canonical commutators between the field operators. The magnetic power 
spectra and the electric power spectra can be computed from the density matrix of Eq. (13. 7p 
by explicitly evaluating the following pair of correlators of field operators: 



,(B) 
'(E) 



Qij '(r,r) 



{I3i{x,T)B,{x + f,T)) = Tr[p5,(f,r)4-(x + f,r)], (3.13) 
{Ei{x,T)Ej{x + f,r)) = TT[pEi{x,T)E,{x + f,T)]. (3.14) 



Since most of the forthcoming formulas are separately valid for magnetic and electric corre- 

(X) 

lators, we shall use the notation Qlj where it is understood that X = B, E. The isotropic 
and time-dependent correlators (r, r) are transverse and can therefore be represented by 
means of the standard fiuid parametrization (see, e.g. [HI SSI 1^): 



(r, r) - 



(r, r) 



or r 



(3.15) 
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where ^^'"^(r, r) and Gi^\r, r) denote, respectively, the transverse and the longitudinal parts 
of the either magnetic or electric correlators; since the electric and magnetic correlators are 
both divergenceless they must also be subjected, to the differential relation appearing in 
Eq. fl3.15p (see also, in a related hydrodynamical context, Ref. [HlllQ]). Using Eq. fl3.14p 
the explicit expressions of the transverse and longitudinal components of the electric and 
magnetic correlators become: 

Ol^\r,T)= J d\nkPx{k,T)q^{kr), g?\r,T) = J d\nk Px{k,T) q^{kr), (3.16) 
where Px{k,T) denotes either the magnetic or the electric power spectrum: 

PB{k, t) = ^|/.(r)|2 + 1), PE(fc, r) = ^\gk{T)\' {2N, + 1). (3.17) 
The functions q^ikr) and qi^{kr) of Eq. (13.161) are given in terms of spherical Bessel functions 



q^{kr) = joikr) 



qhikr) = jo{kr} h J2{kr}. (3.18) 



kr kr 
With the previous notations for the power spectra, Qlj (r, r) can be represented in Fourier 
space as 



1 rd-'k 



Att J k^ 



Px{k,T)P,,{k)e 



-ik-f 



(3.19) 



where Pij{k) is, again, the standard transverse projector. 



3.3 Power spectra from thermal initial conditions 

As discussed in Appendix [Bl Eq. (13.121) can be solved under different approximations. If 
a/a evolves monotonically the magnetic and the electric power spectra of Eq. (13.171) can be 
parametrized as: 



k^ — 

Px{k,T,T^) = Pxix, y, z) = e-2^>l^(x,y)(2Ar, + l) 



(3.20) 



where x{k,Ta), y{k,T,T„) and 2;(r, Tg-) are dimensionless variables: 

x{k,Tcr) = kr^, y{k,T,T^) = / ka{k,T')dT', z{t,t^) = 2tt I a{T')dT'. (3.21) 



The functions M.x{,x, y) with X = B and X = E are computed from the results of appendix 
iBland they are: 

MB{x,y) = \x\H^^\x)\^[ac{y) + i:s{y)f + x\H';^\{x)\'s\y) 



X s(?/)[ac(y) + Es 



H^\x)H^:^\{x) + i7(2)(x)i/?_\(x)] }, (3.22) 



M^ix^y) 



x\H^^\{x)\'[ac{y)-i:s{y)f + x\H^^\x)\\s\y) 



- X s{y)[as{y) - Y.c{y)] 



Hi'\x)H!^\{x) + Hi'\x)Hl'\ 



(1) 



X] 



(3.23) 
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where c{y) = cosh|/, s{y) = sinhy and S = 27ia/k. The time t„ marks the beginning of the 
conducting phase at the end of inflation. In the sudden reheating approximation coincides 
with the end of the inflationary phase. To evaluate explicitly the power spectra we must 
specify N^- The obvious choice iC| 



1 



(3.24) 



where the comoving temperature of the initial thermal background. The power spectra 
f l3.20p contain both the contributions from the vacuum fluctuations and from the thermal 
fluctuations. When A^^ <^ 1 (i.e. k S> /ct) the quantum fluctuations dominate whereas when 
Nk ^ 1 (i.e. k <^ /ct) the thermal fluctuations are the leading source of inhomogeneity in 
the initial conditions. In the limits x <^ 1 and S = 2iTa/k ^ 1, Eq. fl3.20p (for X = B 
and X = E) and Eqs. fl3.22p - fl3.23p lead to the wanted power spectra for the magnetic and 
electric fields: 



PB(k,T) ~ — x|iJ«(^)l'e-'^(^'") coth 
Stt 



PE{k,T) 



k' 



SttS 



2 x\HI% 



x] 



2kT 

2^-2J(k,.) ^^^^ 



r k 1 












2kT. 







(3.25) 
(3.26) 



where J{k, a) = [z{t, To) — y{k, r, Tq-)] which becomes, using Eq. fl3.2ip and the results of 
appendix [B) 



J{k,a) = k 



j:{k,a) - jT.^ik,a) - 1 



dr' 



k' 



47r(j(r' 



-dr'. 



(3.27) 



It is useful to introduce the magnetic diffusivity scale k^^ such that J'{k,a) = k'^/k^ where 
k~'^ = dr' /[inai^r')]. The evaluation of k^- is complicated by the fact that the integral of 
Eq. f]3.27p extends well after r^. This estimate can be made rather accurate by computing 
the transport coefficients of the plasma in different regimes [HUES] (see also [56]). For the 
present purposes, however, such an effort would be forlorn given the minuteness of the ratio 
{k/k„Y the phenomenologically interesting scales. By taking r = Tpq in Eq. (13. 270 the 
result is 

A;\2_ 4.75x10-26 / k 
krr) .foh^oTZTy ZTT VMpc" 



l2hln^Q{z,^ + l) 
where ^mo is the present critical fraction in matter 



(3.28) 



hn is the Hubble rate in units of 



100 km/(sec Mpc) and 2;cq + 1 = cto/c^eq — 3200 is the redshift of matter-radiation equality. 



^The density matrix of Eq. (|3.7I) describes a mixed state obeying the Bose-Einstein statistics. A Bose- 
Einstein multiplicity distribution is not sufficient to infer local thermal equilibrium. In various quantum 
optical situations (see, e.g. [52] pp. 159 and also [l^) chaotic light (i.e. white light obeying as a Bose- 
Einstein distribution for each mode of the radiation field) is generated by a source in which atoms are kept 
at an excitation level higher than that in thermal equilibrium. In this sense the present calculation holds 
also in more general cases when Nk does not have a Bose-Einstein form. 
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Eq. f l3.28p can be also used, with some caveats, to estimate S or fc/S ~ k'^/a. Equations 
(I3.25P and f l3.26p show that the electric power spectrum is suppressed not only inside the 
Hubble radius but also outside by the factor ~ jo^ . The present value of the power 
spectrum can be estimated as 



10'^"-^M o.o ..-J Jf{uXr.M). (3.29) 



Gauss V2.43 x IQ-V V4.15 x IQ-^ 

where P-Qik^r) = a'^{T)P]3{k,T) ancl^ 



f{u, k, T^, k^) = /C(z/) IkT^r'"' coth 



2kr< 



}C{u) = ?^r\u). (3.30) 



The relevant scales of magnetogenesis (say between few tenths of Mpc and 10 Mpc) have 
corresponding wavenumbers which are much larger than the putative effective temperature of 
the thermal background prior to the onset of inflation. The presence of the protoinfiationary 
thermal background implies that 

k 52.54 _(^__^) ^ k \f g,^ y/^f An e ^V4 



A;t (vMpc-V U06.75j I2.43 x 10-V U.Olj ' ^^'^^^ 

where iVmax denotes the maximal number of efolds presently accessible by large-scale obser- 
vations, i.e. 

1 / ^ A / ^0 A , 1 f ^o^Ro 



AT^ax = 62.2 + - In - In — + - In " , . (3.32) 

2 VlO-sy V0.7y 4 V4.15xl0-5/ ^ ^ 

For the numerical estimates it is useful to bear in mind that Tq- = l-L^ = l/[a„H„]] thus /cTq- 
appearing in Eq. fl3.30p is expressible as: 

' 3.22xlO-f^)m-"V.^)"'" (3.33) 



a^H^ ■ VMpc^V VO.Ol/ V2.43 X 10-9. 

As already anticipated in section [1] the value of A^max is close, in the sudden reheating 
approximation, to the minimal number of efolds A^min needed to solve the kinematic problems 
of the standard cosmological model (i.e. A^min — ^max — ^i)- Recalling the fiducial set of 
cosmological parameters determined on the basis of the ACDM paradigm and in the light 
of the WMAP7 data [H iS], Eq. f l332|) gives A^i = 63.6 + 0.25 In e. By definition A^^^ax 
is derived by requiring that the present size of the Hubble radius is all contained in the 
event horizon at the onset of the inflationary phase. The figures given in Eq. fl3.32p are 
for guidance olnly: A^max can be much larger if the reheating is delayed. If right after a 
conventional inflationary phase the Universe expands at a rate which is slower than radiation 
A'max increases. In particular, if, after inflation, the energy density of the plasma is dominated 



^The slope of the magnetic power spectra can be defined, from Eq. ()3.30p . as fc"^ ^ wliere ub = (6 — 2^) 
Note also that, in Eq. jSjQl), /C(5/2) = 1. 
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by a stiff source with sound speed coinciding with the speed of hght we get to the estimate 
Nuiax = 78.3 + (1/3) Ine where it has been assumed that the stiff phase starts after inflation 
and stops right before big-bang nucleosynthesis (see, e.g. [57] and references therein). 

In Eq. (13.311) the value of Qi can be, at most, 0{1). This means that k ^ itl the 
range of scales 0.1 Mpc^^ < k < lOMpc^^ and that the thermal correction is also, at 
most, 0{1), since, in this range, coth [/c/(2fcT)] — C^(l)- The wavelength associated with 
the protoinfiationary thermal background is then always larger than the typical length-scale 
related to the gravitational collapse of the protogalaxy, i.e. (9(Mpc); such a scale is inside 
the horizon prior to matter-radiation equality. Thermal initial conditions are only relevant 
for k < 0.1 Mpc'^ From Eq. flX^ . for u = 5/2 and = A^^^^, ^/P^ik~^ ^ lO'^^ Gauss 

in the range k > 0.1 Mpc while ^JPjs^k, tq) — 10~^^ Gauss for k ~ 10~^Mpc^^ roughly 
corresponding to the present value of the Hubble radius. It is finally appropriate to note that 
for z/ = 1/2 no amplification takes place (since A is constant) and the power spectrum will 
be 0(10^^^) (A;/Mpc~^)^ Gauss: this result is just the value of the magnetic field computed 
from the quantum fluctuations over the protogalactic scale. 



4 Conducting initial conditions 

The most generic protoinfiationary initial conditions are characterized by a globally neutral 
plasma consisting of an equal number of positively and negatively charged species whose 
total energy density and pressure are given by: 

Ptot =P+ + P-+ Pr, Ptot =P++P-+ Pr- (4.1) 

As in the previous section, the bar over a given quantity indicates that the corresponding 
variable is physical (as opposed to comoving). In what follows m± and T± shall denote 
the masses and the temperatures of either the positive or negative charge carriers. As long 
as T± ^ m±, the temperatures and T_ coincide with which is, by definition, the 
temperature of the radiation, i.e. T+ ~ T_ ~ T^. In the opposite case (i.e. for T± < ■m±) 
the evolution of the various temperatures depends on the plasma parameter (?piasma coinciding 
with the inverse of the number of charge carriers contained in the Debye sphere. 

4.1 Thermodynamical considerations 

The different pressures and energy densities appearing in Eq. (14. ip are: 

3 — — 

= m±n± + -n±T±, P± = n±T±, (4.2) 

_ TT^ 4 _ vr^ 4 

Pr = T7 9th T^, Pr = -r7 9th T^. (4.3) 
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From Eqs. f l4.ip . f l4.2p and f l4.3p the first principle of tlie tliermodynamics and tlie adiabatic- 
ity of the evolution imply the following relation: 



d 



Vh (n+m+ + n^mJ) + ^ ( n/T^ + n_T_ 



+ n+T+ + n_T_+pJ(iVH = 0, (4.4) 



where Vh = {^n /'i)H^^a? and is a fiducial volume coinciding, for instance, with the 
volume of the event horizon at the onset of the inflationary stage. Assuming the global 
neutrality of the plasma, the concentrations of the charged species must be equal, i.e. rr+ = 
n_ = TTq and Eq. fl4.4p becomes: 



d[a^{T+ + T_)]+ a-fd{aTr) = 7 



2s 
no' 



(4.5) 



where no = a? uq and s = a^'s are, respectively, the comoving concentration of the charged 
species and the comoving entropy density. Equation (14.50 can be solved under different 
approximations. The physical initial conditions stipulate that, initially, T< T_ ~ Tr with 
the result that the common temperature of the different species scales aqj 



4+7 

T ~ a 2+7 . 



47r^ /n 



7 



45C(3) Krio 



(4.6) 



where rij. = a^riy. and fir = gthT^ C(3)/7r^. If rir ^ no, T scales, approximately, as in the 
opposite case (i.e. n^ -C no) the effective temperature evolves, to first order in I/7, as a~^. 



4.2 Plasma approximation and breaking of duality symmetry 

At finite charge density the gauge-invariant action of the system is 

^tot = Sy-J V^J^^^Y^ d^X + J V^jIT^Y^ d\ + ^charged + ^neutral, (4.7) 

where Sy has been given in Eq. (I2.17P and j^"^^ denote four-currents associated with the 
charge carriers; S'charged = 'S'(+) + S'(_) and S'neutrai denote, respectively, the actions of the 
charged and neutral species. By minimizing the action (14. 7p the evolution equations can be 
obtained and then presented in terms of the normal modes introduced in Eq. f l2.19p : 



-Lv ■ {^/XE) = 47rg(n+ - n_), v^V ■ 

1 Id 
-^V X {y/\B)=ATiq{n+v+-n^v^) + — — (VaE), 



(4.8) 
(4.9) 

(4.10) 



^Recall that C(3) = 1.202. 
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wherj^ q = e/\/X. The normal modes of Eq. fl2.19p have been used in section Eland in 
appendix |A] in the source-free case. It is then natural to employ the same variables also at 
finite charge density. 

The source terms in Eqs. (14. 8 p and (14. 9 p break the duality symmetry discussed in Eqs. 
(I2.20p and (I2.2ip . The main difference between quantum and classical initial conditions re- 
sides in the breaking of the duality symmetry. In the case of the quantum initial conditions 
the duality symmetry is broken at the end of the inflationary phase. Conversely for classical 
(or conducting) initial conditions the duality symmetry is broken already during the protoin- 
fiationary phase. Even if Eqs. (14. 8 p and (I4.9p -( 14T0|) hold in general, for direct comparison 
with the results of section [3] the attention will be focussed on the case of homogeneous \/X. 

Assuming an approximate kinetic equilibrium between the charged species, the first of 
Eq. (14. 8 P implies — = 47rg(?2+ — ?2_) with n±{(f)) = noexp [^q(j)/T±\ where denotes the 
Coulomb potential; expanding n±((f)) for \q(f)/T±\ < 1, we get — = 0. The explicit 
form of the comoving Debye screening length is 



T+ T_ IT 
— )■ 



, v^. (4.11) 

V 47rg^no(-/+ + J-j V 87re^no 

The second expression of Eq. (14.111) defining £d holds in the limit of approximate thermal 
equilibrium between the charged species, i.e. T+ ~ T_ — )■ T. In the presence of a background 
of charged species the correct expansion parameter is not the gauge coupling itself but rather 
the inverse of the total number of particles present inside the Debye sphere: 



1 4 T 

^plasma = ND = -nnoe%, = \ h 5 ; (4.12) 

5'piasma plays the role of coupling constant of the plasma and it measures the degree to 
which collective effects dominate over single particle behaviour. In other words the plasma 
parameter is small when many particles interact at the same time. For the Debye shielding 
to occur the number of particles in a Debye sphere must be large. But this means (/piasma ^ 1 
and the latter condition defines the conventional plasma approximation. 

The validity of the plasma approximation ensures that the system can have a charge 
density, a shielded electric field but very limited interaction of single particles. It is easy to 
check that 7 appearing in Eq. (14. 6 p can be written in terms of l/sfpiasma showing, as antici- 
pated, that the corrections to the adiabatic scaling of the effective temperature, proportional 
to 1/7, are indeed small in the plasma approximation. If Weyl invariance is unbroken the 
plasma parameter is also Weyl invariant. Since the plasma parameter does not explicitly 
depend on the masses of the charge carriers the following chain of proportionality relations 
holds: 

(^plasma OC (riQ il)-^ CC (^q T^)'^/^; (4.13) 



^"According to the rescaling established m the previous section, we should have that q — > (7/(47r) in 
Eqs. (|4.7|) and (|4.9|) . We prefer, however, to redefine e and keep the Gaussian units where the 47r appear 
throughout since these conventions are the standard ones for the analysis of weakly coupled plasmas. 
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but tiq is insensitive to the expansion as it is evident by passing from comoving to physical 
quantities: (noT^)-^^^ = {n^T )-^/^. The Weyl invariance of the plasma parameter implies 
that the initial value of (y'piasma will be preserved by the dynamical evolution. 

If Weyl invariance is broken the validity of the perturbative expansion and the increase 
of the Debye shielding scale can be used as criteria for enforcing the consistency of the 
perturbative approximation and the consequent screening of all electric fields throughout 
quasi-de Sitter stage. According to Eqs. (14.111) and (I4.12p £d must increase with -\/A since 
only in this case the electric fields will be progressively screened as the plasma evolves 
during the protoinfiationary phase. But this is indeed the same requirement necessary for 
the decrease of fifpiasma- In the opposite case (i.e. when \fX decreases) the electric fields are 
not properly screened so that the Debye scale at the end of inflation is actually smaller than 
the value it had at the onset of inflation. 

The validity of the plasma approximation implies a preferred evolution for the Debye 
shielding length. It is therefore possible to trade the evolution of for the evolution of £d- 
On a phenomenological ground, it seems reasonable to adopt, as working guess, the following 
expressior0 

hm ; \ = J — — , hm ; \ 

for a ^ Or, £d — ^ constant, the plasma is again Weyl invariant (at least as long as the charge 
carriers are relativistic) but duality is still broken at finite conductivity. In the case of quasi- 
de Sitter dynamics, Eq. f l4.14p reproduces, by construction, the evolution of a/A already 
discussed in section [3] in the complementary case of quantum and thermal initial conditions. 
This coincidence facilitates the comparison between the two sets of initial conditions. 




(4.14) 



4.3 Lorentzian plasmas and conductivity 

Even if -\/A does not evolve in time, the Weyl symmetry is broken because of the mass of the 
charge carriers. In section |2] it has been shown that Ohmic currents are Weyl invariant in 
the relativistic limit where masses and chemical potentials are negligible in comparison with 
the temperature. The contribution of the charged species to the transport coefficients can 
be computed within a kinetic model of the plasma. Working under the hypothesis that the 
collisions between the particles of the same charge can be neglected, it is consistent to assume 
that V± ^ r+ and F-t r_, as it happens for Lorentzian plasmas [HH |59]; furthermore, at 
high temperatures T±> H for interactions mediated by massless gauge bosons. 

Note that denotes the conventional normahzation of the scale factor at the reheating epoch. The 
evolution of the Debye shielding increases during the protoinfiationary phase and also during the inflationary 
phase. The intermediate exponential suppression at reheating may be absent but has been included to 
account for a possible breaking of the plasma approximation in the transition regime. 
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The Vlasov-Landau equation for the positively and negatively charged species is given 

by: 

^ + ■ ^'^^ + ^ ■ ^^^^^ = l^j.on' ^'-''^ 
where tit and p± are, respectively, the comoving three-velocities and the comoving three- 
momenta defined as: 

^ = ±q[E + v^xBl v^= . ^ (4.16) 

If P| ^ m^a^ the whole system is invariant under Weyl rescaling and the evolution of the 
geometry does not affect directly the evolution of the charged species. In the opposite limit 
the two relations appearing in Eq. fl4.16p can be combined and the equations for the velocity 
field become: 

v'+Hv+='^[E + v+xB]+ aT^^{v_ - v^) + \P^aT+r{vr - t/+), (4.17) 
p+a p_ 3 p+ 

^_ + 'Hv^ = -^[E + v^xB]+ aT^^{v+ - v^) + l-^ar_,(tv - v^), (4.18) 
p_a p_i_ 3 p_ 

where is the velocity of the radiation fiuid whose dynamics is not relevant for the present 
ends; r±j. are the rates of momentum exchange between charged and neutral species while 
r± are the rates of momentum exchange among the charged species themselves. 

To deduce the expression of the conductivity Eqs. fl4.17p and f l4.18p must be combined 
to obtain the evolution of the center of mass velocity and the equation for the total current: 

v = ■ , J = qno{v^ -v_). (4.19) 

m+ + m_ 

The result of the mentioned combination is: 

^' + = w/''^/ V + t^^^r-{Vr - V), (4.20) 

a*p(l + m_/m+) 3p 

J' + iH + aT± + -P — - ]J =-^iE + vx B + 

V iriQ m^J 47r V g no 

4gp^r^_ ^ 

H — [V-Vr). 4.21 

3m_ 

where p = + and Wp is the plasma frequency. To simplify the expressions a hierarchy 
in the masses of the charge carriers can be assumed (for instance m+ > m_ = m). In Eq. 
fl4.2ip the terms J' and 1-LJ are comparable in magnitude and are both smaller than T± and 
T±r. Thus the explicit form of the Ohm law becomes: 

J = a(E^,.B^^-i2ll\ (4.22) 

V qriQ noq / Anal ± 
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where the contribution of the neutral species has been neglecteco. For a Lorentzian plasma 
[58| |59] , the conductivity can then be written as: 

a(a, 7) = — p£i2)_, lim T(a, 7) oc a'^^^^^^^ = constant. (4.23) 
?V^ + ^ 

Consider first the case when the comoving Debye length does not evolve (i.e. £d constant or, 
in equivalent terms, z/ = 1/2 in Eq. (14.141) ). In the limit 7^1 the comoving temperature 
T = Ta is approximately constant and this is verified provided fifpiasma ^1- In the limit 
T ^ ma Eq. fl4.23p implies cr ~ T/g^ (as it happens in the case of a relativistic plasma); in 
the opposite limit, a ^ T /q^^T / {ma). 

When ^D evolves as in Eq. f l4.14p . the breaking of the Weyl invariance produced by 
the finite mass of the charge carriers interferes with the evolution of the Debye shielding. 
After expressing Eq. (I4.23p in terms of the number of inflationary efolds N = In (a/a^,), the 
evolution of the conductivity becomes: 

JN+Ni){2u-1) 

^(N) = a, . ; (4.24) 



VT+e 

note that Ni accounts for the initial number of efolds during the protoinflationary phase 
while Nc is given by: 

1 1 / m \ 

= -0.253 + - In e - - In (7th - In (^—j . (4.25) 

Following the notations of Eq. fl4.14p . Eq. f l4.24p holds before reheating, i.e. in the regime 
where initial conditions are set. To leading order in (^piasma the critical number of efolds 
A'^c depends on the temperature reached during the protoinflationary phase and it can be 
estimated by recalling that at the onset of the inflationary phase the total energy density 
does not exceed the contribution of the protoinflationary energy density stored, for instance, 
in relativistic species. With this logic, from Eq. f l4.25p . an upper bound on can be derived 
and it can be expressed as 

iV, ~ 36.78 - 0.25 In (^th/100) + 0.5 In (^/lO"^) - In (m/GeV). (4.26) 

Given the value of Nc there are four physically distinct situations. If ~ Atot ^ Ac + A'min 
and u = 1/2, Weyl invariance is broken before the onset of the last (9(63) efolds of inflationary 
expansion; this means, in practice, that the sources do not contribute to the initial conditions 
which are accurately fixed by quantum mechanics. If A^min < A^ < Ac + A^min and z/ = 1/2 
the conductivity will be constant for the first A"c efolds and then it will be exponentially 
suppressed as e^^--^--)/^ ^ _ jy^.j (jf _ jy^.^ + jyj_ If ^ ^ 3/4 



^^In what follows the Hall term and the thermoelectric term will be considered to be of higher order and 
they will be anyway irrelevant for the present ends. 
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the conductivity is practically constant for N ^ Nc, Finally if u ^ 3/4 and z/ 7^ 1/2 the 
evolution depends on the specific value of z/. 

The goal here is not to endorse (or predict) a specific duration of the inflationary phase 
but just to convey the message that when A^min < N < Nc + N^ij^ the last C(63) efolds 
of inflationary expansion may start when the conductivity did not undergo a substantial 
suppression. For instance, when the mass range of the lightest charge carrier is 0{GeV) and 
if ~ (9(A^jnin) the conductivity is still almost constant (9(30) efolds prior to the end of 
inflation. In this class of physical situations the normalization of the electric and magnetic 
fields does not follow from the quantum mechanical initial conditions but rather from the 
conducting initial conditions. 



4.4 Plasma initial conditions 

Equations (I4.23P and (I4.24p imply that for N < Nc the conductivity can be assumed to 
be roughly constant only if z/ = 1/2. When u ^ 1/2 Eqs. (USD, (USD and (HTTU]) must 
be numerically integrated to determine the evolution of the electric and magnetic fields. 
Using the normalized scale factor as evolution parameter Eqs. (14. Sp . (14. 9 p and (I4.10p can be 
reduced to the following pair of vector equationj^ 

f = + ,4.27) 

da aT-i aV. oH 

f = ^8-^, (4.28) 
da aH aTi 

where a = a/a^,; a^, denotes conventionally the value of the scale factor when the inflationary 
phase starts and T-L is given by: 




n = n.M-] + - • (4-29) 



According to Eq. (I4.29P the protoinflationary phase is dominated by a perfect fluid with 
barotropic index w; the inflationary phase is simply modeled by an effective cosmological 
constant. More complicated models of protoinflationary evolution can be formulated by 
introducing scalar field sources together with perfect fluids. A class of analytic solutions 
of the Friedmann equations with this property is reported in appendix O For the present 
purposes, the details of the inflationary sources are immaterial. This has been explicitly 
checked by using the solution illustrated in Eqs. ( IC.ip and ( ]C.2p . 



Following the previous discussions and, in particular, Eq. (I4.14p Eqs. (I4.27P and (I4.28P 
can be integrated conductivity evolves as 

a{a) = ^= ^ ^ , C = e (4.30) 
y'l + a/a, C 



^■^ Equations ()4.27|) and (|4.28p are supplemented by the requirement that W-E = W- B = Qa.s iniphed in 
the case of a globally neutral electric plasma. 
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Figure 1: The result of the numerical integration for the case w = 1/3. 



( depends on the specific value of A'c. By arbitrarily assuming that a is constant for N < Nc 
we can expect that the ratio between the electric and the magnetic power spectra approxi- 
mately obeys the scaling law: 



valid in the limit a ^ a,,,. The numerical results confirm this guess. In Figs. [T]and|2]the 
system (14.271) and (I4.28P has been integrated for two illustrative values of the barotropic 
index, i.e. w = 1/3 and w = 1. On the vertical axis, in both figures, the common logarithm 
of the ratio between the electric and the magnetic power spectra is reported; on the horizontal 
axis we have the common logarithm of the normalized scale factor. The value of if* is given 
in Planck units while a^, is essentially the normalization of the conductivity in units of H^:. 
In both figures the dashed, dotted and dot-dashed lines correspond to three different values 
of K = k/H^, which is the wavenumber in units of H^. The full (thin) lines correspond to 
the analytical approximation of Eq. (I4.3ip valid for but plotted also at earlier times 

just to guide the eye. 

In the limit where all the fields appearing in Eqs. (I4.27P and (I4.28P are solenoidal (i.e. 
V-.E = V- B = V- J = 0) the displacement current can be neglected for very high 
conductivity and therefore the appropriate initial conditions for the electromagnetic fields 
at Tx are simply given by 



as assumed, without explicit numerical proof, in [H]. The power spectra subjected to the 




(4.31) 





(4.32) 
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o.=10'', H. =10"*', C= 10"", V =5/2 




Figure 2: The result of the numerical integration for the case w = 1. 



conducting initial conditions f l4.32p and subsequently amplifield according to Eqs. fl4.14p 
and dMZD-dOHD arS 



PE{k,T,T,)=PB{k,T,)\V{z,x) 



(4.33) 



where x = kr^, z = —kr and 7] = k/Aira; the expressions 

in fx" ( 1 



V{x,z) 



U{x, z) 



4 V z Ix 



Pl^\z)Pl^\x) - Pl^\x)Pl'\z) 



Hi'\x)P^'\z) - Hi'\x)P^'\^ 



ITT [Z 



2u — TjX 



Hl'\x)Hi'\z) - Hl'\x)Hi; 



2)/ 



+ X 



(4.34) 



are given in terms of the Hankel functions of first and second kind (i.e. H^\z) and H^\z)) 
as well as in terms of the following combinations: 

Pl'\z) = 2uHl'\z) - zHIU^), Pl'\z) = 2uHl'\z) - zHIU^). (4.35) 



By setting r = —t^ the expressions of Eq. (14.331) reproduce the conducting initial conditions 
given in Eq. (I4.3ip . Expanding V{x,z) and U{x,z) in the limit a; -C 1 the following 
expressions can be obtained: 



V{x,z) 



yjVi(x,z) + ,/^V2{x,z) + o(^x^/^^ 



(4.36) 



""^^The dimensionless variables x and z introduced hereunder and appearing in Eqs. (|4.33p . ()4.34p and 
(j4.35p must not be confused with the dimensionless variables of section |3l Since the two sets of variables 
never appear simultaneously, potential confusions are avoided. 
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z) 



X 



(4.37) 



where the following functions have been introduced: 



Vi{x,z) 

V2{X,Z) 

Ui{x,z) 

U2(X,Z) 



) — 5— f 



xJu{z) 



x{rix -2)T{u - 1) +4r/r(z/) 



TT J_^(^)[32z/(z/ + 1) - lQr]{u + l)x - 8{u + l)x^ + Ar]x^ + x^ 
A/xr(z/ + 2) sin (ttz/) 
2^^-=^ J,_i(z)[4r/r(z/) + r(z/ - l)a;(r7a; - 2)], 
2-5-^^ TT Ji_^(z)[8r(z/ + 2)((x + 2r7)x - 4z/) - V{u + + 4r/)] 

v/ir(z/ + l)r(z/ + 2)sin (ttz/) ' 



(4.38) 



where Ju{z) is the ordinary Bessel function. For \tx\ ^ |Te| (and x < 1, z < 1) the 
electric power spectra are suppressed throughout the whole stage of inflationary expansion. 
Note that the rate of the suppression is larger than in the quantum and thermal cases 
examined in section |3l From Eqs. f l4.36p . f l4.37p and (14.381) . at the end of the inflationary 
phase, PB{k,Te,T^) ~ PB{k,r^){ae/a^)^''-^ and P^{k,T^,T,j.) ~ r]PB{k,T^) {ae / a^)^'"^" . Since 
7] = k/{4TTa) ^ 1 the suppression of the electric fields is always much larger than in the case 
of vacuum initial conditions. 




Figure 3: The bounds of Eq. (I4.40p are illustrated for Nc = 35 and A^max = 65. 

The amplitude of the magnetic power spectrum during the protoinflationary phase can 
be maximized by only asking the compatibility with the closure bounds. In the latter case 
the power spectrum of the magnetic field will be given by 



_ 10-60.15 '°^-^^(.-l/2)(^-^.)-2(^-iV_) .4 

Gauss llO-V \hJ ' ^ ^ 
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where S = y PB{k,Ti)/{HfMp) < 1 measures the fraction of energy density stored in the 
magnetic field at ti, Hr accounts for the possibihty of a delayed radiation-dominated phase 
between the end of inflation and the onset of big-bang nucleosynthesis. The exponent ai 
depends on the expansion rate between the end of the inflationary phase and the onset of 
the standard (i.e. post-inflationary) radiation-dominated epoch. Equation fl4.39p has several 
interesting limits. In the case u = 1/2 there is no amplification due to the evolution of 
the gauge coupling and therefore the upper bound on the magnetic field intensity is around 
0(10"^^) Gauss in the sudden reheating approximation where Hr ~ Hg. It is interesting that 
this figure coincides with what has been obtained in Eq. fl3.29p for u = 1/2, k = 0.1 Mpc~^ 
and ^ A^max- This result coincides with the magnetic field one would obtain from the 
protoinflationary initial conditions in the case of standard thermal history and minimal 
number of efolds. If z/ = 5/2 and iV = iVi ~ (9(65) and ~ (9(35) the maximal magnetic 
field turns out to be 10~^^ Gauss (in the sudden reheating approximation) which can become 
of the order of O{10^^^) Gauss for a stiff post-inflationary phase extending down to the 
nucleosynthesis scale. 

The classical and the quantum results can be more quantitatively compared by imposing, 
for instance, the following reasonable hierarchy of inequalities: 

^^/Mk^\ ^ /^^(^x ^ / y^Pg(A:,ro) x ^ ^^^^^^ 

V Gauss / quantum V GaUSS / classical V GaUSS / maximal 

The quantum and classical contributions both have to be smaller than some maximal value 
which is computed by assuming that the magnetic power spectrum is of the order of the 
energy density of the inflaton at the end of inflation. The results of this comparison are 
reported in Figs. |3] and H] for different values of A^^c and A^max- On the vertical axis the 
value of the normalized rate is reported and a function of the total number of efolds, i.e. 
5'(A^) = J-'/'H which coincides with (z/ — 1/2) during the quasi-de Sitter stage of expansion. 
With the full line the bound stemming from the maximal value of the magnetic power 
spectrum is indicated. For values of >S'(A^) and A^ within the shaded region the hierarchy 
expressed by Eq. f l4.40p is satisfied. Below the dashed line, in both plots, the first inequality 
of Eq. fl4.40p is inverted and the quantum contribution is larger than the classical one. 
Figures [3] and m are purely illustrative and they show that both classical and quantum initial 
conditions lead to sizable magnetic fields. The differences between the two sets of initial 
conditions must come from a more thorough analysis of the spectral properties especially for 
typical length-scales larger than the Mpc. 

5 Monopole plasma 

The duality symmetry can be used to deduce what happens if instead of an electric cur- 
rent we have a monopole current. In a globally neutral monopole plasma the value of the 
total magnetic charge vanishes i.e. = n_ where now n+ and n_ denote the monopole 
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Figure 4: The bounds of Eq. (I4.40p are illustrated for Nc = 10 and A^max = 85. 

concentrations with opposite charges. The masses of the positively and negatively charged 
monopoles are assumed to be different (i.e. m_|_ ^ mJ). There is no possible confusion with 
the notations used for the electric plasma since we are not going to consider here dyons 
but simply a globally neutral plasma of magnetic monopoles. In the past the possibility 
of a dyonic plasma in magnetohydrodynamics has been discussed in [601 El]; the situation 
considered here is, by construction, even simpler since only magnetic charges will be consid- 
ered. Even if this is strictly not necessary, the monopoles are assumed to be non-relativistic 
both during the protoinflationary phase and during the inflationary phase. The transport 
coefficients can be computed by exploiting the assumption of Lorentzian plasma [58] in full 
analogy with what has been discussed in section H] but with the crucial difference that the 
role of the magnetic and of the electric degrees of freedom is somehow interchanged. The 
relevant system of equations is given by: 



where e denotes the monopole charge. The explicit form of Eq. (15. ip becomes then: 




V^(AZn = 0, 



<4l =±e^n^n< 




(5.2) 




(5.3) 



(5.4) 



(5.5) 
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where Q = note that v± denote the velocities of monopoles with opposite charges. From 
Eq. 05. 2 p the evolution of the comoving three-momenta can be derived and it is given by 



dP± 
dr 



±q 



B-v±x E 



R 



v± 



± 



(5.6) 



In the relativistic limit ^ im?a?, v± = P±/\P\ ■ In the case P^ <^ m^a^ the evolution 
equations of the velocity will be manifestly not invariant under Weyl rescaling: 

4 Pr 



vV + 'Hv+ 



^[B-v+xE] + af^^iv. - v+) + -£^ar. 
P+a P- 3 p+ 



(5.7) 



— [B -v^x E]+ ar^ — (v+ - v) + -^aT^Jvr - V-)- (5.8) 
p_a p+ 3 p_ 



Defining the center of mass velocity of the monopole system and neglecting the momen- 
tum exchange between the monopoles and the radiation background the following pair of 
equations can be obtained: 



JjnX E 
«^Pm 



— # — # — * 

Jra = (^m{B - X E) , 



(5.9) 



where Jm = rim^v^ ~ '^-) is the total monopole current and cr^ is the magnetic conductivity. 
The main equations of the electric and magnetic fields in the presence of the monopole 
current are given by 



y/Xdr -dJVXEj 



A 

Xd. 



B 

71 



V'^E 



V^B 



-47rV X J„ 



An 



drUx.l 



(5.10) 
(5.11) 



We can assume, as before, that y/X is homogeneous. The situation is dual to the case of 
the electric charges: the magnetic fields are screened while the electric flux is conserved at 
finite magnetic conductivity. The induced magnetic field will then be -B ^ —V x E /{ATram)- 
It is interesting to speculate that an initial monopole plasma could screen the magnetic 
components but not the electric ones which could be instead amplified if \/X decreases 
(to keep the magnetic field screened); the electric components may be converted back to 
magnetic fields at reheating. The analysis of these themes is beyond the aims of this section 
whose only purpose is to illustrate the role of the duality symmetry in connection with the 
initial conditions of inflationary magnetogenesis. 



6 Concluding remarks 

Which are the correct initial conditions for protoinfiationary magnetogenesis? Are they clas- 
sical? Are they thermal? Are they quantum mechanical? Are plasma sources immaterial 
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when setting the initial conditions of inflationary magnetogenesis? Are the obtainable mag- 
netic fields phenomenologically relevant? These are some of the main questions addressed in 
the present investigation whose motivation stems from the observation that vacuum initial 
conditions for the evolution of large-scale gauge fields are customarily imposed regardless of 
the number of inflationary efolds and in spite of the possible presence of protoinflationary 
remnants. It will be recorded that conventional inflation is known to be geodesically incom- 
plete in the past and while it is true that the contribution of the sources to the evolution of 
the geometry is likely to be exponentially suppressed with the number of efolds, the same 
conclusion does not apply to the evolution of large-scale gauge fields. 

The symmetries of inflationary magnetogenesis suggest that the most generic initial con- 
ditions are neither quantum nor thermal but rather conducting: if charged particles are 
not flne-tuned to vanish, the simplest situation compatible with the phenomenological con- 
straints is to contemplate the case of a globally neutral protoinflationary plasma. Conven- 
tional inflation takes place when the gravitational coupling is strong. As inflation proceeds, 
the temperature of the protoinflationary plasma decreases at a rate which depends on the 
smallness of the plasma parameter reflecting the largeness of the Debye shielding scale which 
must increase if magnetic fields are to be amplified and electric fields are screened. In this 
situation the electric fields are suppressed in comparison with the magnetic fields not only 
asymptotically in the future but also at the onset of the protoinfiationary phase. 

The Weyl and the duality symmetries determine the initial conditions of infiationary 
magnetogensis since, as usual, the symmetries of the problem are reflected in the symmetries 
of the solutions. Quantum initial conditions do not break explicitly the duality symmetry 
but may break the Weyl symmetry if the gauge kinetic term is coupled to a spectator 
fleld. Classical initial conditions may not break the Weyl symmetry but break explicitly 
the duality symmetry since they imply, in the simplest case, the presence of Ohmic currents 
without a corresponding magnetic source. The present considerations demonstrate explicitly 
that both conducting and quantum mechanical initial conditions lead to phenomenologically 
relevant large-scale magnetic fields so that possible distinctions between the two sets of initial 
conditions must rely on subtle differences in the spectral properties for typical length-scales 
much larger than the one of the protogalactic collapse. 
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A Duality- invariant Hamiltonian 



The duality symmetry leaves the Hamiltonian invariant but when the Hamiltonian depends 
explicitly on time (as in the present case) every canonical transformation changes its form. 
Which is the correct Hamiltonian to use for the quantization and for the evolution of the 
electromagnetic fluctuations? The punchline of the following appendix is that there indeed 
exist different Hamiltonians, all differing by a canonical transformations, but only one class 
of Hamiltonians is explicitly invariant under duality and this will be the preferred one, for 
the present purposes. 

In the Coulomb gang , and in fiat space-time, the duality properties have been discussed 
in [21 [3] where it has been noted that even if it is possible to keep both Yq and the longitudinal 
part of the Abelian vector potential (i.e. Il), the transverse variables are decoupled from 
the longitudinal ones and from the gauge contributions. Without any gauge fixing, the 
extremization of the action with respect to Fq implies that Y^^' = VYq (see |21 E])- Thus the 
action of Eq. f l2.17p can then be written, in the Coulomb gauge, as 

Sy = J dTLyir), Ly{t) = J d^x Cy{x,t), (A.l) 

with 

£y(x,r) = + J'V - ■ y' - dS- d'y]. (A.2) 

The canonical momentum conjugate to y can be obtained from Eq. (lA.2p and it coincides, up 
to a sign, with the canonical electric field defined in Eq. fl2.19p . i.e. vf = y' — J^y = —E. The 
canonical Hamiltonian is simply obtained from the Lagrangian density Hy{t) = J d^x [vf ■ 
y — Cy{x, t)] and its explicit form is: 



HY{r) = ^Jd^x 



n'^ + 2TTf-y + diy-d'y 



(A.3) 



The Fourier mode expansion for the canonical fields 



^(^' = / d'k vf.^r) e-^^-', y{x, r) = / d^k U-) ^-^'^ (A.4) 



can be inserted into Eq. (1A.3P and the resulting expression is: 



Hy{t) = lJd'k 



^k-^-k + J^[^k- v-k + ^-k-yk] + k% ■ y 



(A.5) 



The Hamiltonian (lA.SP is invariant under the transformation provided, at 

the same time the electric variables are appropriately rotated into the magnetic ones and 
vice-versa, i.e. 

k^-^Vk^ Vk^l^k^ ^-k^-ky.k^ y-k^l^-k, (A.6) 



TT 



^^The Coulomb gauge condition is preserved under Weyl rescaling of the time-dependent metric; the 
Lorentz gauge condition does not have the same property; see the discussion after Eq. (|2.5p . 
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where k = \k\. The Lagrangian given in Eq. f lA.2p is not invariant under the duahty 
transformation of Eq. (1A.6P and this happens since the duahty transformation maps a 
tensor into a pseudo-tensor. The Hamilton equations derived from Eq. (]A.5|) become: 



re, 



(A.7) 



Under the transformation of Eq. ( 1A.6|) the two equations of Eq. ( 1A.7|) are transformed 
one into the other and vice-versa. Since the Hamihonian (]A.5|) is time-dependent its form 
can be changed by a canonical transformation. The Hamiltonian (1A.5I1 is invariant under 
the transformation (1A.6P exactly because the canonical momenta are, up to a sign, the 
canonical electric fields. If the Hamiltonian is different, the explicit duality invariance is not 
guaranteed. Consider, indeed, the following functional of the old fields (i.e. y), of the new 
momenta (i.e. H) and of the conformal time: 



n[y,n,T] 



(T X 



y-Ii 



7 



-y ■ y 



(A. 



The old momenta are related to the new ones as vr = H — J^y while the new Hamiltonian 
becomes Hy{t) ^ iJ^^'^V) = Hy{t) + d^n, i.e. 



(new) 
Y ' 



r) 



(Tx 



If the calculations are carried on in terms of H 



(new) 
Y ' 



(A.9) 



(and not, as preferable, in terms of 



Hy{t)) the spectra of H will not be related by duality transformations to the spectra of -B = 
V X y simply because H does not coincide with the canonical electric field E. It is therefore 
mathematically convenient and physically justified to use Hy{t) rather than ify"™''(r). In 
the literature this point is never mentioned but it seems essential when computing the 
magnetic and the electric power spectra with either quantum or thermal normalization. 



B Mode functions: explicit expressions 

The solution for the mode functions can be discussed with different techniques depending 
on the monotonicity properties of A. If A has a monotonic dependence on the conformal 
time coordinate, the evolution equations for /^(r) and 5'A;(r) can be solved exactly. Suppose, 
for instance, that \/X = V^i(— t/ti)^/^^*^. Thus, the rate of variation of is given by 
J-" = [1/2 — v)/t and Eq. fl3.12p can be solved in terms of Hankel functions: 

fk{r) = v^if«(-fcr), Af = y|e-(^+V2)/2^ (g.^) 
9kir) = -MJ^V^H^\{~kT). (B.2) 
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The dual of fk, Qk and J-' are, respectively, fk — )■ Qk/k, Qk — > —kfk and J-' — )■ — J-". Duality acts 
non-trivially on the solutions given in Eqs. (IB.ip and (lB.2p . Indeed, if J-" — )• J-" = (1/2 — /i)/r, 
= —J-" provided z/ = 1 — /x. Duality can then be used, ultimately, to relate the spectra of 
the electric and magnetic fields in the ideal (and somehow unphysical) situation where the 
finite conductivity effects are absent both at the end and at the beginning of inflation. 

Whenever the evolution of a/A is not monotonic or not solvable analytically in a closed 
form, Eq. f l3.12p can be separately solved in the limits k/J^ > 1 and k/J-" < 1. When 
k/J" > 1 



fk{r) 



9k{r) 



-jfc(r-rin) 



(B.3) 



valid for r < Tex where, by definition, Tex is J^{tcx) = k. In the limit k/J^<l the solution of 
Eqs. f l3.12p can be formally expanded in powers of the spatial gradients: 

oo oo 

fk{r) = E k''hAr), 9k{r) = E k'^gkA^), (B.4) 
£=0 e=o 

with the result that fkj and gkj obey the following hierarchy of coupled equations: 

fk,i = 9kj + J'fk,e, i>0, (B.5) 
9k,o = -J^ 9k,o, ^ = 0, (B.6) 
9k,i = -fk,e-i - 9k,e, ^ > 1- (B.7) 

Equations (1B.5I) - (]B.7|) are analytically solvable order by order and for a generic form of the 
pump field. The lowest order solution with the boundary conditions of Eq. (IB. 31) is simply 

1 



fk{r) 



9k[r) 



'A(r) 
Apx 



Ik 



X 



A(r 



cx ^-ik{n 



in) ■T(r T 



dr' 



(B.8) 



(B.9) 



where Aex = A (rex). When the Universe reheats, the conductivity of the plasma breaks the 
duality symmetry and the relevant equations are given by: 

9k = -k^fk - 47r(T^fc, = gk. (B.IO) 

Defining with and gk the solutions of Eqs. fl3.12p the solutions for r > can be obtained 
rather easily by direct matching with the result that 



fk{T) 



kfki^a 



ka{k, a 
+ ^fc(^^)sinh[?/(r,r^)] 

g-2(r,To-) 



a{k, a) cosh [?/(r, r^)] + S(/c, a) sinh [y{r, r^)] 



fB.ir 



gk[r) 



gk[^a) 



a{k, a 

- ^7fc(T-<T)sinh[?/(r,r^)] 



a{k, cr) cosh [?/(r, r^)] - S(A;, a) sinh [t/(r, r^)] 



(B.12) 
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where 

a{k, r) = ^j^\r)-l, S(r) = (B.13) 

Note that S(r) measures the ratio of the conductivity over the wavenumber. In Eqs. fIB.lip 
and (IB.12P y(r, t„) and z(r, t„) depend on the evolution of the conductivity and are defined 
in Eq. ([MI])- From Eqs. ( IRTT]) and (IRT2|) \fk{r)\^ and |fi'fc(r)p are easily obtained and 
they can be studied in the limit k <^ 27r(T (i.e. S ^ 1) where the conductivity dominates. 
The opposite case (i.e. when k ^ 27rcr) is obtainable from Eqs. (IB. lip and (1B.12P by 
appreciating that, in this limit, a ^ if3, y ^ iy and z ^ iz where (3 = a/1 — S^. In the 
latter situation, the expressions of |/a;(t)P and [(/^(r)^ are different since a, y and z become 
all complex quantities. It is finally appropriate to remark that the evolution of the canonical 
operators is such that [yi,T^j] — > for r > t„. Because of the presence of the conductivity 
the Wronskian is driven to zero. This means that out of the two solutions of the system 
only one survives, i.e. the one related to the magnetic part. The vanishing of the Wronskian 
signals the transition to the classical dynamics where the magnetic field operators become 
Gaussian random fields. 



C Protoinflationary evolution: analytic example 

An explicit solution describing the protoinflationary dynamics is [62] : 

3{w + 1) 



a{t) = 



sinh(/3if^t) 



1/^ 



cp{t) = ipo±\I^Mp^Jl-n, In 



tanh 



2 

V 2 



(C.l) 
(C.2) 



where w is the barotropic index characterizing the protoinflationary fluid while a{t) and (p{t) 
are the scale factor and the inflaton expressed in cosmic time; the parameter fi* measures 
the fluid fraction of the protoinflationary energy density 



ft. 



3Hm 



2 ' 



Ptot(^) = P*(^^ 



(C.3) 



where ptot denotes the protoinflationary energy density. The potential for ip is given by 



■{l-w)H^Ml{l-Q,)smh^ 



1(3 ((^ - <fo) 



2 {i-n,)Mp 



(C.4) 



Equations (IC.ip and (lC.2p satisfy the solution of Eq. (13. ip . For f3H.t < 1 the solution is 
decelerated and from Eq. (IC.ip we have a{t) ~ a.{(3H.tY/^ where Hi = 2/[3(w + l)ti]. 
In the opposite limit (i.e. l3H.t ^ 1) the solution is accelerated with H{t) ^ since 
H{t) = if,/[tanh iPHj)] and H = -PH^/[smh.^ (PHJ)]. 
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